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The thermodynamic properties of nuclei are studied in a mean field model using a Skryme in- 
teraction. Properties of two component systems are investigated over the complete range of proton 
fraction from a system of pure neutrons to a system of only protons. Besides volume, symmetry, 
and Coulomb effects we also include momentum or velocity dependent forces. Applications of the 
results developed are then given which include nuclear mechanical and chemical instability and an 
associated liquid/gas phase transition in two component systems. The velocity dependence leads to 
further changes in the coexistence curve and nuclear mechanical and chemical instability curves. 
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I. INTRODUCTION 



One primary goal of medium energy nuclear collisions is a detailed study of the thermodynamic properties of 
strongly interacting nuclear matter [JQ- An important feature of these properties is the existence of liquid-gas phase 
00 ' transition. Properties of the nuclear force (long range attraction and short range repulsion) parallel those of a van 
der Waals system [3] which qualitatively describes a liquid-gas phase transition in atomic systems. The liquid-gas 
phase transition in nuclei is the first phase transition seen in a strongly interacting system. Relativistic heavy ion 
i ' collisions are being used to explore a second phase transition from hadronic matter made of mesons and baryons 
to a quark-gluon phase. Important differences exist between the nuclear interaction and interaction between atoms. 
Because nuclei are made of neutrons and protons, the phase transition is in a two component or binary system where 
symmetry energy effects and Coulomb effects play an important role. Moreover, the nuclear force has a velocity 
dependence. The presence of symmetry energy and a Coulomb interaction effects and also a velocity dependence in 
the nuclear interaction makes the nuclear case a unique and interesting binary system within the general scope of such 
systems. Examples of other two component systems are binary alloys and liquid 3 He. For 3 He the two components 
are spin-up and spin down fluids. The phase structure in such two component systems has some important features. 
In nuclear systems, isospin fractionation [l], 0, 0, [H, S 0] is an example where the monomer gas phase has a large 
neutron to proton ratio. Ref.[7] is the most recent reference to isospin fractionation and contains further references to 
, it. Both the symmetry energy and Coulomb energy play an important role in this phenomena of isospin fractionation. 
' Nucleons carry spin but very little research has been done in understanding the role of spin in the liquid gas phase 
\ structure. However, the crust of neutron stars has features associated with a superfluid phase. 

An early study of the nuclear liquid-gas phase transition [3J] treated the system as a one component system of 
nucleons. This study was then extended to two components using a Skyrme interaction Q. A relativistic mean field 
model was also developed in Ref.@ where the role of the symmetry energy was studied in detail. The addition of 
the Coulomb energy (lO . [Tl| resulted in asymmetries which changed the mechanical and chemical instability regions 
and binodal surface in pressure P, temperature T and proton fraction y associated with phase coexistence. For 
one component systems a phase diagram is the more familiar binodal curve of pressure versus density or volume 



determined by a Maxwell construction. Some other studies of one and two component phase transitions can be found 
in Refs. [H Ha, Q QJ, QJ, [H, [H HI ■ The present work is an extension of our research reported in Refs. 0, El • 
An extended Skyrme interaction is now used in our present study and, for example, includes effects associated with 
a velocity dependence in the nuclear interaction. Here we will use a simplified form of the velocity dependence. 
In particular we will use an effective mass approximation for it which includes a density dependent behavior. Our 
primary goal is to see what qualitative effects a momentum dependence has when superimposed upon an interaction 
model that does not include them. A momentum dependence study was also given in Ref.[22| using a more refined 
dependence for it. Our results differ from that of Ref.[22| since we also include Coulomb and surface effects. Coulomb 
effects lead to an asymmetric behavior in proton fraction [ToL of various quantities. In the absence of Coulomb 
forces a symmetry exists around proton fraction y = 1/2. The velocity dependent force modifies nuclear saturation 
properties and the symmetry energy. Some recent extended studies of symmetry energy can be found in Refs.f23l.l24l]. 
The results given below show modifications in chemical and mechanical instability curves arising from an inclusion 
of a density dependent effective mass. The velocity dependence has a larger effect on the proton rich instability and 
coexistence features compared to the neutron rich curves. A detailed discussion is given below in Section IIIII and in 
associated figures. The study of two component nuclear systems with arbitrary neutron/proton ratios will be useful 
for future RIB (Rare Isotope Beam Facility) experiments and in astrophysical studies such as in neutron stars. 
Our paper is divided as follows. The next section discusses the thermodynamic properties of nuclei. It is divided 
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into two subsections. General results based on a mean field approach are presented in IIA. Specific results based on 
a Skyrme force for the potential terms, and low and high temperature kinetic energy behavior appear in IIB. This 
subsection also contains the effects of a velocity dependent interaction and related effective mass results. Then in 
Sect lllll we apply the results of subsection IIA and IIB to the specific issues of: A) mechanical and chemical instability 
of nuclei and B) the liquid-gas coexistence curve. Results are presented in 9 figures which are discussed. Finally, in 
Sect. IV a summary and conclusions are given. 



II. THEMODYNAIC PROPERTIES OF NUCLEI IN A MEAN FIELD DESCRIPTION 



A. General Results 



In this section we present results for the thermodynamic properties of nuclear matter which are extended from the 
results of Ref. [101] to include a velocity or momentum dependent interaction. The matter is a two component system 
of protons and neutrons in equilibrium at some temperature T. We first develop expressions for the total energy E as 
a function of the density p q of each component q and temperature T. The behavior of the energy functional with p q 
and T can be used to obtain the behavior of the pressure P and chemical potential p q for each component of type q. 
These quantities will also be functions of p q and T . They can then be used to study, for example, a phase transition 
in the nuclear system. 

To begin, we use the fact that at a given temperature T = 1//3, the proton and neutron constituents are distributed 
in phase space according to the Wigner function / as 

f(r,p) =^/ g (f ? ,p), f q (r,p) = ^f q (r,p) = ^ 1 (1) 

9 

The spin degneracy factor 7 = 2 and e q and p q are the single particle energy and the chemical potential of particle of 
type q. Then the particle density p and nucleon number A are given by the following equations: 

pif) = E pi ' pi ^ = / d3 pf* ^ p) > ( 2 ) 

q J 
A = Y^N q = J d 3 rp(r), N q = J d 3 rp q (r) = J ' d 3 r J d 3 pf q {f,p) (3) 

Defining r(r) as 

r(f) =Y,r q {r), T q (f) = J d 3 p^f q (f,p) (4) 

the total energy E is given by 

E = J d 3 r£(f) = J d 3 r J d 3 p^-f(r,p) + J d 3 r J d 3 pU(f,p) = J d 3 r [8 K {r] + U (f)} (5) 

The potential energy density is U(f), while the £k(t) = j^ T ('' f ) i s the kinetic energy density. The single particle 
energy e q is given by 

SE 58(f) p 2 6U p 2 

The u q — jj- is the single particle potential of particle q which may in general be momentum dependent. The chemical 
potential p q is given by e q at an effective Fermi momentum p = pp q defined by the following equation: 

2 

N = eq\ P=PFq = -^ + Uq(r,PF q ) (7) 

In order to study a phase transition we need information about the behavior of the pressure when the system 
is in equilibrium. The general expression for the pressure can be defined dynamically from the total momentum 

conservation law, ^ [J d 3 r J d 3 p pf] = — J d 3 rV r - 11= 0, using the Vlasov equation as developed in Ref. plj : 

^ + (V p e 9 ) • (V r / 9 ) - (V r e 9 ) ■ (V p / 9 ) = (8) 
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A more general expression is obtained from the hydrodynamic consideration of TDHF in phase space as given in 
Ref.[26| which reads; 



v r - n 



d_ 

dt 



d 3 pp^2f q (r,p) 



E/ d V 



dfq 
dt 



£ / d 3 ppV r ■ [(VpeJ/J + E/ d3 PP- (y r e q )f q 



(9) 



where p — p/p is a unit vector in the direction of p. Using (V r e q )f q = V r (e 9 / g ) — e 9 V r / 9 = V r (e g / g ) — V r £, the 
dynamical pressure tensor 11^ is given by 

% = E / d3 pp* ( v p £ ?) /« + / rf3 ?E e ?/« - £ 
= E / d3 ^ v £ 

9 ^ 

= E / d3 ^ 



5S_ 



fq + Sij 



5U 



E/^ 



6S_ 



fa-£ 



fq + h 



E/^ 



SfqJ Jq 



u 



(10) 



Our previous study [101 ] focused on a momentum independent potential which gave the following simpler results for 
the pressure tensor: 



% = E/ dV-^L 



E 



su 

^Pq 



Pq-U 



= /^Efv«+^E 



The diagonal element of Tlij is the pressure P — Tlij which simplifies to 

p = n« = E / rf3 A? + E - 17 = Pa ' + E^« - ^ = Pk + Pp 



(11) 



(12) 



The Pr- = J d 3 p^f = \&k is the kinetic part of the pressure P, while the interaction potential part is Pp = 
£<? u qPq ~ U = p 2 S ^ U J p P ^ ■ At temperature T = the pressure P is related to the derivative of the energy per particle 



E I A with particle number fixed as 



P = IL 



dV 



dp 



(13) 



This result applies to a single component system. Below we will give results at non zero temperature for a multi 
component system. We first proceed with a discussion of the role of the momentum dependence and effective mass. 

As mentioned, our study is based on a qualitative study of the role of a momentum dependent interaction and 
we therefore use a simplifying approximation. Specifically, we use an effective mass with a density dependence and 
this approximation greatly simplifies our analysis in two component asymmetric and finite nuclear systems. We still 
include both Coulomb and surface effects since realistic nuclear systems have such terms which are important in their 
description and stability properties. More refined stiudies will be developed in future work. When the momentum 

2 

dependent part is of the form A(p p , p n )fr/i then it can be incorporated into the Hamiltonian as an effective mass 
term. In Ref . [22j ■ the momentum dependence is obtained from 



P P l + (p-j7)7A2 



(14) 



We use an effective mass m* / m approach (2f| for Eq. (fT4| which can further be approximated by expanding the factor 
1/(1 + (p — p 1 ) 2 1 A 2 ) to first order in 1 — (p — p 1 ) 2 / K 2 . Specifically, we write the effective mass behavior of m/m* as 



m 

m* 



a , \ 2m 

i + A(p)tt 

n 



MP)Tq(?) 



2m 



2m 



r q {r) 



(15) 
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Moreover, we have U(r) — U(p) + A(p)T(r) = U(p) + ^ q A q (p)r q (r) with r q (r) of Eq.Q. A momentum dependent 
single particle potential u q (f,p) is given by 



"'^ " ^(r,p)- ffp, + S Pq + A ^tf- S Pq 



p 2 SU(f) i ^ / 



The is related to u q (f,p) through the result 



P-q 



P F q 

2m 



+ u q (f,p Fq ) = l+A q {p) 



2m\ p 2 Fq SU{p) 8A{p)r{r) 



2m Sp q 



1 + A q {p) 
a 



2m\ p\ q SU(f) _ Pp q 8U{r) 



2m Sp q 2m* Sp q 



(16) 



(17) 



Also 



and 



J d 3 pu q f q 



V p u q (r,p) 



5U(r) 
~6f q ~ 



d 3 p 

SU(r) 
Sp q (f) 



611(7) 



mp) = 5 M^ Pq{ f) + 5 ^^ Pq (f) + A q (p)r q (r) 



p q {r) + A q (p)T q {r) 



Here p and r are treated as independent variables. Then the pressure tensor fly is given by 



(18) 



(19) 



E 

L q 



f 5U{r) 



Pq + A(p)r(r) - U(r) 



r -* ,^ 6(U(p)/p) M(p)r(r) 



SPq 



SPq 



p(rlPq(rl S -^P^ + A(p)r(r) 



q "1 q 

and the pressure P or diagonal element 11^ = P is 



&Pq 



^ = n, = E(i + ^)|?)/^ + E 



p 9 + A(p)T(r) -C/(r) 



p 9 - = P* K + P P 



(20) 



(21) 



TheP£ = E g ^ g with 



* f = (l + ^)?f)/^/.-/^/.-|^ 



(22) 



is the kinetic pressure with an effective mass correction term, and the second equality is for an isotropic momentun 
distribution. The potential part of the pressure Pp is given by 



8(U(P)/P) , SA(p)r(r) 

PPq 7 1" Pq- 



5Pq 



5p q 



E 



PPq t — +A q (p)T q (r) 
°Pq 



p2 W(p)/p) , p 5A(p)r(f) = p2 S(U(r)/p) 



6p 



Sp 



6p 



+ A{p)r { r)=/-^M + r K - £K 
Sp 



(23) 



5 



The = ^(f) and 
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Also in obtaining this result we use the fact that 



«-Esr« ( "-s( 1+ ^) 



2m \ ?i 
2m 



r 9 (r). 



Y.pi — - — - = /' 



Sp 



(24) 



(25) 



which can be shown by looking at the derivative of B(p)C{p p )D{p n ) = B(p)C(px p )D(px n ). Here the variation p must 
be done after replacing p q by px q . 

Other thermodynamic variables, such as S, ft, F, G are given in Ref.[lGj. The entropy S follows, from the 
distribution f q of Eq. (JTJ , 



and 



S 



S q = -jL I d*p 



/ 9 ln/ g + (l-/,)ln(l-/,) 



(26) 



^ f f q -Pp q jd pf q 

In equilibrium, from Eqs. fTU)) and (j2"7) 

TS = £ + P - 2_j ^iPl = £ K + Pk - ^(Pq ~ Ug)pq 



(27) 



Pig 

2m 



Pq, 



(28) 



The last eqaulity of Eq. (f28|) is the result of using Eq. ([7|) . For momentum dependent potential the entropy is now 

Pq 



TS = £ + P-^p jqPq = g* K + P * K -^L q - 5 J^\ 



£ 



K 



Pf 9 



2m 



2m? 



(29) 



where use has been made of Eqs. pTl) and (|2"Tj) to obtain this result. General thermodynamic relations also determine 
the entropy, pressure and chemical potential [10( . 



B. Thermodynamic Properties of Nuclear Matter based on a Skyrme Interaction 



We now use a Skyrme interaction to develop expressions for the potential U. Once the potential energy U in Eq.© 
is known, then questions related to mechanical and chemical instability and the possibility of a phase transition of the 
system can be studied using Eqs.([T]) - (|2"9"|) . The potential energy U determines e q and p, q and the potential energy 
part of E and P. Then for fixed T and N q , the Wigner function / and pp q are determined and thus the kinetic terms 
of E, fi q , and P. Using these results, the entropy S can be determined. For a nuclear system of proton (p p ) and 
neutron (p n ), this gives the local potential energy density as 



U(Pq) 



t 



2 V 2 J H 



1 



to 
2 



l r 

~4 



U 1 



2 



ti (l 



Xo 

2 



£<* + 5( 1+ t)^ 



PT 



h 

12 

1 

2 



'■2] ^PqTq 



+c P l3 p 2 p + c a p* 



(30) 
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TABLE I: Skyrme parameters used here are in MeV and fm units For ti and to., the SkM parameter values are used. 



to x 


h x 3 


a 


-1089.0 -1/6 


17480.4 -1/2 


1 




Momentum-Dep. 


Momentum-Indep. 


h 


251.11 





XI 


-1/2 


-1/2 


t 2 


-150.66 





X2 


-1/2 


-1/2 


Effective mass m* jm 


0.895626 


1 


Binding energy E B /A 


13.1057 


15.54447 


Fermi energy Ef 


31.8018 


34.2101 


Saturation density po 


0.1283 


0.143145 


Symmetry energy Sv 


23.4791 


24.39379 


Compresibility k 


307.780 


361.9045 



Here Cp 13 = ^-e 2 R 2 and Csp* 1 = 47rJ ^ / °'(p) = I^Ifl£lp 2 /3 wnen we approximate the Coulomb and surface effects as 

coming from a finite uniform sphere of radius R = r^A 1 ^ with total charge Z (Uc = f^^f^) 03- The values for the 
force parameters used here are given in Table [J We define an effective mass m* as 



m 

m* 



2m ( 1 



xi 



U 1 + — )+t2 1 + — 



■r-2 



ti - + Xi 



Then the momentum dependent potential term becomes 



MP) 



1 



2m* 2m 4 



2m 



-1 + 1 



(l ( 1 + a )+i2(l + ?) 



X-L 



*2 ( g + X 2 



X2 



2m f 1 r 



si - to 



Pq 



Pq 



X2 



(31) 



P<i 



For a symmetric nucleus, N — Z , p q = p/2, and thus 

„a+2 



U(p) = h a p 2 + ^t 3 p a + 2 + 1 (tt + t 2 ) pr + C(Ppl + C sP v 



(32) 



(33) 



This potential enegy determines the interaction dependent terms of £, P, e q , and p, q which depend on densities without 
an explicit T dependence. 



For a momentum dependent potential energy as in Eq. (|30p , e q — p q = (p 2 — p 2 Fq ) / (2m* q ) where the effective mass 



2 2 

i* is independent of the momentum independent part of potential and the Wigner funcion of Eq. |T]) becomes 

1 1 



fq(r,p) 



of + 1 e 0(p 2 -p 2 Fq )/(2m* q ) + x 



(34) 



Thus we can evaluate the kinetic terms in £, P, and p q which are functions of T and PFq- Defining the Fermi integral 
F a (r]), with effective mass m*, 



P Pq 



1 + e *-V q 

5U(f, 



-dx 



AttTi 2 



a+l 



2p 2a+1 dp 



o 1 + e 



p 2 /2m*-77 9 ' 



SPq 



>lq 

X q = J2nh'/m q T 



Pp Fg /(2m*) = p 2 Fq /{2m* q T) = lnz„ 



we can write, for f(r,p) = f(r,p), 



Pq = J d 3 P fq(r,p) = ^ j dV 



I 27 

/3(p 2 - P ^)/(2m;) - " = \ 7^ Fl /2(^. 



(35) 

(36) 
(37) 

(38) 



7 



~Fq 



PFq _ ft (tor 



2m* 2m* \ 7 



v3 P 



-Pq 



2/3 



d V~^fq{r,v) - h3 



7 / 



d 6 p- 



1 



H 2 e /3(P 2 -P 2 F ,)/(2m*) + 1 



6k 2 
2m \ 7 



2m* 



-Pq 



2/3 



8 7 V^A- 5 F 3/2 ( % ) = -_/^4A- 3 F 3/2 (ry g ) = ^£ Kq = ^£* Ka , 



r 3 

2^ = 2^ 



d 3 P^fq(?,P) h3 

m* 1 27 



47 ft 2 v/tt 



7 / * P 



2m e /3(? 2 -P^)/(2m;) + 1 



c* _ 
C Kq — 



3 



2m; v 2' K? 
— — : — \ ^3/2 W 



m /3 
2 



T q = ^-PrTa = / d 3 P ^-fq{f,p) 



7 / .3„ F 



K 3 F 3/2 (Vq) 



2m^ e /3(p 2 -p| g )/(2™j) + x 



1 2 7 ,_ 3 



(39) 



(40) 



(41) 



(42) 



Here tFq is the chemical potential at absolute zero or Fermi energy and pFq is the effective Fermi momentum at T 
(which is related to density p q through Eq. (|55|) ). The particle number N q = J d 3 rp(r) determines the effective Fermi 
momentum pp q (r) or rj q at T, in terms of density p q (r), 



Vq(pq,T) = (3 [fig 



_ su(f) \ p^ ! (VI 



SPq 



2m* V2 



V2 7 



A; 



Pq 



(43) 



For multi(two)-component systems with potential energy given by Eq. (|30|) . with a given p q (or PFq) and T, the 
thermodynamic properties are as follows. The chemical potential is given by 



fJ- q (Pq,T) = Ti] q (p q ,T) + 



SU(r 



T Vq ( Pq , T) + t a (l + ^-)p+^(l + ^)(a + 2)p 



a+l 



12 V 2 



to o + x o ) Pi + 



*3 ( 1 



12 V 2 



+ x 3 (a-l)2p a p q 



h (I 



12 V 2 



*1 1 + TT + *3 1 + TT 



J' 2 



1 

T 

4 



+ x 3 ap 



+ x 3 2ap a - L p' 



1 



*1 - + Xx - t 2 - + X 2 



+Cf3p f) - 1 p 2 p + 2Cp p p 6 q>p + vCsp*'- 1 . 
The equation of state has a behavior determined by 

p/ rp\ v-2 2 S(U(P)/P) SA(p) T (r) 

P(Pq> T ) = 2^^ £ Kq{P^ T ) + P ~ +P 



(44) 



E 
E 



<5p 

£Kq(Pq> T ) ~ £Kq{Pq,T) 



dp 

S(U(r)/p) 



Sp 



o £ Kq(Pq^ T ) " £Kq(Pq,T) 



to 



x 0\ 2 *3 



12 



■''3 



E^"fl^ + ^)(a + lKE 



12 V 2 



+C(/3 + l)/^ + C s (^lK. 

The energy density is 

ffo.T) = J2 £ Kq(Pq,T) + U(f) =^2^^) + U( P ) 



(45) 
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+ !( 1 + f)p 2 -!G-o)Ep? + &(i + f)^-SG-^"E^ 

\ / q \ / q 

+Cp f} P 2 p + C* s p" (46) 
and the entropy can be obtained from 

TS{ Pq ,T) = ^^/feTj-^-^K^^rl-T^^^TK. (47) 
q q q q q 

Once we evaluate 1*1/2(77) and -£3/2(77), or more directly i] — F^ix) and F^^iv): we can evaluate various thermody- 
namic quantities in terms of p q and T. 

For low temperature and high density limit, A 3 p large, i.e., when the average de Broglie thermal wavelength A is 
larg er than the average interparticle separation p _1//3 , we can use a nearly degenerate (Fermi gas) approximations 
[27f for F 1 / 2 to obtain 



Vq(Pq,T) 



B(u m W\-(j& q _p-l(^ A 3 fl 

^ ^7 J - ^ - F v 2 V« 



= /3e 



7T 

1 , 

12 \ e 



T 



2m; V 7 
2 7 



2/3 



2/3 _ n " l q 

Pq 



J_ 

ifi A V6tt 2 
3 



4/3 



A rv- 2 / 3 + ... 



(48) 



£Kq(Pq, T ) = 7T^ X q ^3/2 (ty?) = T^Kq = ^Pq e *Fq 



3ft 2 /(hr 2 ^ 2/3 
10m* \ 7 



1 



t>7T 

~L2~ 



5tt 2 to* 2 



3h 



4 V6vr 2 / 



4/3 



TV/3 



(49) 



2m* _ 3 2m; , 
h 2 Kq ~ 5 ft 2 Pq€F * 



1 



5tt 2 / T 



3 /6tt; 

7 



2 \ 2/3 



12 

5tt 2 ( T' 



-Fq 



3 /e^ 2/3 

5 ^ 7 



12 



5?r 2 m* 2 



3?i 4 V67r 2 y 



(50) 



In the other limit where X q p is small, we have a nearly non-degenerate Fermi gas (classical ideal gas) and the resulting 
equations are given by an ideal gas in leading order with higher order corrections [13] as 



Vq(Pq,T) 

£-*Kq\Pq-> T) 

T q(Pq, T ) 



Pq 



SU(r 
Sp q 



= In 



T^Kq ~ ^PqT 

2m *q P * _2m q *3 

2~^Kq — 



PqXq 

7 

1 PqK 



1 Pq^ 



h z 2 



2 5/2 



PqT 



7 



2V2 7 

1 2 

8 ~ W 2 



In 



Pq 



1 / PqXq 

2V2 \ 7 



1 P 9 A 



2 5/2 7 



V 7 

1 2 

8 ~ W 2 



P9A 3 



(51) 
(52) 

(53) 



For a nuclear system with protons and neutrons with the interaction given by Eq. (j3"0|) , the non-degenerate Fermi 
gas limit of Eqs. (|51"j) . (j52")l and ([53")) leads to the following set of equations. The chemical potential has a behavior 
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determined by 



tl 1 + 



■I'l 



2m! 



Pg 



3 2m| 
2 7i 2 



Pg 



25/2^,^9 



+(o(l + ? ) p+ i|( 1 + |) (£t + 2K «_i|(' + , 3 )^« 



1 



+ #0 Pg + 



*3 /I 



12 V 2 



+ x 3 )(a-l)2p a p q 



t 3 (I 



12 V 2 



a; 3 2ap' 



Q-l 2 



The equation of state has a form given by 



5 T 
22V5 ■ 



Ag 



*3 



3 x - 



g 
m 



Pq 



2V2 \ 7 



Pg 



n a+2 



2 \2 



*3 f 1 



12 \ 2 



+CC8 + V^ + C.^-lip". 



The energy density is 



£{p,y,T) = 
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The effective mass m* and thus X q arc, in general, isospin dependent [251 ]. However we will consider an isospin 
independent effective mass here for simplicity in this present study. For the case of m* — m* with A g = A (such as 
the case of x\ = Xi = — 1/2), these equations become: 
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Here, for the proton density (p p ) and neutron density (p n ), we defined, the isoscalar density p, isovector density P3, 
proton fraction y and related quantities by 



p = Pp + pn, P3 = Pp~ Pn = (2j/-l)p, y = p p /p, 

Pp = Up + Ps) = yp, p« = |(p - P3) = (1 - y)p, 
E 9 P 2 ^ = \(P 2 + Pi) = [1+(2 r 1)2 V = [l + 22/(2/ - i)]p 2 , 

^ - ip(p 2 + 3pd = ^ +3(2 r l)2 y = [1 + 3^ - i)]p 3 



(62) 



e 9 p 3 , 

The ± in p q stands + for q =proton and — for neutron. 

At fixed T and P, only one of either p or y is the independent variable. Thus observables such as P, £ / p, S / p may 

have a discontinuity in T or y when (^§§^j or (f^) diverges. We can study the behavior of thermodynamic 



quantities at a fixed P using dP = from Eq. ([55|) . 
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This equation gives ue{p) where both dP/dy 
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The yE(jA is indepenent of p for momentum independent Skyrme interaction with x% — —1/2 and (3 = as considered 
in Ref . flQ. [TTj | . The X3 term and the density dependent effective mass for a momentum dependent Skyrm force 
introduce a small p-dependence in yE- Eq. (l59|) shows that, for p-dependent yE, the -P(p) curve for different values 
of y at fixed T may cross at some p. Moreover, the minimum pressure for a given T and p (i.e., (dP/dy) p .T = 0) 
occurs at y = yfi(p) ^ 0.5 due to Coulomb effect. These results were not seen in Ref.@. At yE, the pressure of 
the coexistence curve is minimum and the liquid and gas phases have the same proton fraction yE- The condition 
dP/dy = determines the equal fraction point yE- 



III. APPLICATIONS TO NUCLEAR MECHANICAL AND CHEMICAL INSTABILITY AND THE 

LIQUID-GAS PHASE TRANSITION 

A. Mechanical and Chemical Instability 

We now use the results to discuss features of the instability of nuclei, both mechanical and chemical, and the 
liquid-gas phase transition. The region of mechanical instability is determined by the condition 



dP 
dp 



= (65) 

y,T 



FigJT] shows the behavior of the pressure P(p,y,T) as a function of p for several values of the proton fraction y. All 
curves are at T = 10 MeV. The range of y is from y = 0, or pure neutron matter, to y = 1, or pure proton systems. 
The point y = 1/2 corresponds to symmetric systems. Without a Coulomb interaction results would be symmetric 
about the point y = 1/2 which would also be the point of equal concentration in a liquid/gas phase coexistence. 
Including a Coulomb interaction shift the equal concentraion point to a proton fraction of y = ys(p) ~ 0.415 with 
a momentum dependence included in the interaction and to y = yE = 0.41057 without a momentum dependence. 
The y = yE curve for both the momentum dependent and independent cases has the lowest pressure versus density 
dependence, i.e., the lowest P for a given value of p at a given T. A higher or lower y raises the pressure at a given 
density. Both a momentum dependent Skyrme interaction and a momentum independent Skyrme interaction results 
are shown for several values of y and they are distinguished by the thickness of the lines as described in the figure 
caption. The momentum dependence increases the pressure in the range shown (thick lines compared to thin lines) 
and introduces the density dependence of Ve{p) given by Eq. (|64p . The mechanical instability densities for each y 
curve at T — 10 MeV are the points where the P(p,y,T) curve has zero slope, dPjdp\ y ^ = 0. The total region of 
mechanical instability is obtained by a similar calculation of P(p, y, T) at different T. For a one component system 
or a symmetric system the mechanical instability region is a curve somewhat similar to an inverted parabola with its 
peak at the critical point. Allowing for systems with different values of y gives a two dimensional boundary surface 
for the mechanical instability region. The intersection of the surface with different y planes gives the one dimensional 
boundary curve or line of mechanical instability for each corresponding value of y. 

Fig[2]shows the proton fraction y versus the density p for different fixed values of the pressure at a fixed temperature 
of T = 10 MeV. The loops and curves are determined by solving P(p,y,T) = P for the values of P listed in the 
figure caption and at the temperature T = 10 MeV. Fig. 2 is obtained from Fig.l by drawing a horizontal line and 
looking at the points where the horizontal line intersects the set of P(p, y, T) curves. This intersection can be at one, 
two or three points. Besides the inner most closed loop (P = 0) shown in Figf^l a vertical line exists at p = for 
P = for all y = ~ 1. Similarly, for the second inner closed loop at P = 0.015 MeV/fm 3 , a nearly parallel vertical 
line is present at very low density. The right most point on each curve and the left most point on a closed loop 
with dp/dy\p t r — are at the point of equal concentration yE- Also shown are two thin lines for j/s(p). The dashed 
thin line is at yE = 0.41057 and is horizontal or density independent and corresponds to the momentum independent 
interaction. The solid thin line is nearly horizontal with a slight density dependence and has j/b(p) = 0.4106 ~ 0.4214 
for p = ~ 0.15 fm -3 . Horizontal turning points on each curve occur at dyjdp\p^ = 0. For each T, there is a curve 
P(p)\ y ,T with an inflection point for a particular y which we call yi- At the pressure P = P(p, yi,T), the closed loop 
in Figf2]just breaks at the point of y = yE on the left low density side and creates two new horizontal turning points 
with dy/dp = 0. Fig. 2 also shows the result that a momentum independent force has closed loops outside those of 
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Density, p [fm 3 ] 



FIG. 1: Pressure P(p) versus p at T = 10 MeV for various proton fraction y. The upper solid curve is for y — 0, dashed curve 
for y — 0.2, dash-dotted curve for y = 0.5, dotted curve for y — 0.8, and dash-dot-dot-dotted curve for y = 1. The thick curves 
are for the momentum dependent Skyrme force and the thin curves are for the momentum independent Skyrme force. The 
lower thick solid curve is for y = ys(p) of Eq. (|64[l (j/b = 0.4106 ~ 0.4214 for p — ~ 0.15 fm -3 ) with momentum dependent 
Skyrme force and the lower thin solid curve is for y = JAe = 0.41057 with momentum independent Skyrme force. 



a momentum dependent force and open curves to the right of those of a momentum dependent force with the same 
pressure P. 

The region of chemical instability (spinodal in fJ*(y)\p,T) is determined by the condition 



dy 







P,T 



for each component q = p or n. These conditions for either protons or neutrons give the same relation since 

ydfi p + (1 - y)dfi n = -dP 
P 



(66) 



(67) 



13 



c 
o 

o 
p 



c 
o 

2 

Q. 




Density, p [fm 3 ] 



FIG. 2: Proton fraction j/(p) for P = 0, 0.015, 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5 MeV/fm 3 from inside to outside at T = 10 MeV. 
The solid curves are for the momentum dependent Skyme force and the dashed curves are for the momentum independent 
Skyrme force. The thin straight lines are JAe(p) for the corresponding force. 



This general condition will be used later in our discussion of results given in various figures. The result is also useful 
for checking numerical results. The chemical instability condition can be rewritten in terms of derivatives of the 
chemical potential and pressure with respect to the density variable p and proton fraction y. Namely, the chemical 
instability condition can be obtained from the following relation [llT | 



dP 


dp q 




dP 


dp q 




dp 


y,T d V 


P,T 


dy 


P ,T d P 





(68) 



The expressions developed for the proton and neutron chemical potentials are functions of the variables (p, y,T). The 
equation of state P(p, y, T) can then be used to find their behaviors in terms of (y, P, T) or (p, P, T). The behaviors 
with y of the proton chemical potential p p (p, P.T) — > p p (y) and neutron chemical potential p n {p, P.T) — > p n (y) at 
various values of the pressure P and at a fixed temperature T = 10 MeV are shown in FigJ3J The chemical instability 
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FIG. 3: Chemical potential p p (y) and fi n (y) for P = 0.015, 0.05, 0.1, 0.2, and 0.5 from top to bottom curve for protons (solid 
curve) and from bottom to top curve for neutrons (dashed curve) at T — 10 MeV. Thick curves are for momentum dependent 
Skyrme force and the thin curves for momentum independent Skyrme force. 

region boundaries are determined by the points where the slope of each chemical potential with respect to y is zero. 
Further discussion of the chemical spinodal line is given in the next subsection. The behaviors of the proton chemical 
potential p, p (p,P,T) — > p, p (p) and neutron chemical potential p, n (p,P,T) — > (i n (p) with density p at various fixed 
values of the pressure P and at a fixed temperature T = 10 MeV are shown in FigJU Fig|5] and Fig0] show some 
similarities in the behavior of the plotted quantities, i.e., inner closed loops at low pressure, to outer curves that 
almost form closed loops with increasing pressure, to open curves with further increases in pressures. 
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FIG. 4: Same as Fig[3]but for fip(p) and Hn(p) versus p. 

B. Liquid-Gas Phase Transition and the Coexistene Curve 

For a one component system the coexistence curve is a line obtained by the familiar Maxwell construction as 
already noted. For a two component system the coexistence region is a surface obtained as follows. The condition 
for coexistence between the two phases requires the proton chemical potentials to be the same in two phases and, 
similarly, the neutron chemical potentials must be the same in the two phases at a given pressure and temperature. 
Note that the proton fraction need not be the same in each of the two phases. In fact, the li quid phase should be a 
more symmetric system than the gas phase because of the symmetry potential as seen in Refs.jlOl. Illj. Figs show 
features of the coexistence curves together with the mechanical and chemical instability curves. 

The condition of phase coexistence corresponds to a rectangular box geometrical construction in the chemical 
potential plots of Fig [3] or of FigHl Namely, the chemical potential equality condition p, p (yi, P, T) = fJ, p (y2, P, T) 
and n n (yi,P,T) = ^i n (y 2 , P,T) leads to a rectangular box in Fig[3]with vertical sides connecting the n P {y\,P,T) to 
the p, n (yi, P,T) for side 1 and the fj, p (y2, P,T) to the p, n (y2, P,T) for side 2. The horizontal sides are the chemical 
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FIG. 5: Pressure P versus proton fraction y for coexistence loop (thick curves) at T — 10 MeV. The solid curve is for the 
momentum dependent Skyrme force and the dashed curve is for momentum independent Skyrme force. The thin curves are 
the chemical instability boundary curves for each case of Skyrme interaction respectively. For both momentum dependent 
and momentum independent cases the maximum of the chemical instability loop and the coexistence loop occur at the same 
point where the curves are tangent to each other as discussed in the text. The point of equal concentration is j/_e ~ 0.415 for 
momentum dependent case and ys = 0.41057 for momentum independent case. 



potential equality conditions at y\ and ?/2 for neutrons and for protons. The rectangular box shrinks in its horizontal 
direction in fi q -y plots as the point of equal concentration, where the liquid and gas phases have the same proton 
fraction, is approached (the lowest point of the coexistence curve in Figj5]). 

FigfS] shows various features of the coexistence region in pressure versus proton fraction. The coexistence region 
are the dark thicker solid line for a momentum dependent force and the dark thicker dashed line for a momentum 
independent force. Also shown are associated chemical instability regions as a thinner solid line, and thinner dashed 
line. The calculations are done at a temperature of 10 MeV. For a two component system, the coexistence and 
instability regions are two dimensional surfaces in pressure, temperature and proton fraction as mentioned above. 
The pressure-proton fraction behavior shown is a consequence of cutting these surfaces with a constant temperature 
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plane. The result at T = 10 MeV are the loops shown. Other temperatures can be obtained in a similar fashion. 
For a momentum independent force the chemical instability region basically lies inside the coexistence curve and 
peaks at the top of the coexistence loop, the critical points. The condition dP/dy\T = with d 2 P/dy 2 \T < gives 
a critical point on the coexistence curve and the condition dy/dP\x — gives the point with maximal asymmetry at 
the left and right most points of coexistence curve. The proton rich y > yE and neutron rich y < yE loops are very 
asymmetric because of the Coulomb interaction. The inclusion of velocity or momentum dependent interactions leads 
to further modification of the coexistence curve and chemical instability curves. This modification is easily seen in 
the figure by comparing the dashed momentum independent curves with the solid momentum dependent case. The 
figure shows that the momentum dependent interaction that was used has a larger effect on the asymmetric proton 
rich loop (jj > yE) significantly reducing its maximum pressure. The maximum of the neutron rich loop (y < ys) 
remains somewhat unchanged with a small increase. Another effect is to shift the two loops inward toward the equal 
concentration point yE- A third effect is to shift the lowest pressure point, which occures at the equal concentration 
yE, upward with the value of yE nearly unchanged. Finally, it should be noted that the peaks of the coexistence and 
chemical instability curves are at the same point where the curves are tangent to each other. We see no indication of 
a truncation effect in our model where the coexistence curve intersects the chemical instability curve before reaching 
the peak critical point. A truncation effect gives a limiting pressure (below the maximum pressure of the chemical 
instability curve) above which a liquid-gas phase transition cannot take place (22J. 

Fig[S] shows plots in y versus p of phase coexistence curves, instability boundary loops for both chemical and 
mechanical instability, and features of d/i q /dp\ yi T = for proton and neutrons. The thin curves are for a momentum 
independent interaction and the thick curves are for a momentum dependent interaction. The calculations are done 
at a fixed temperature of 10 MeV. Some features common to both cases are as follows. The mechanical and chemical 
instability boundary curves are closed loops with the mechanical loop (dashed line) inside the chemical instability 
loop (solid line). These two loops touch at yE, the dash-dot-dot-dotted line. The yE(p) increases slightly with p for 
a momentum dependent interaction while it is constant (horizontal) for a momentum independent interaction. The 
p dependence of yE(p) come from the p dependence of the effective mass and also from the X3 term as can be seen 
in Eq. ([M|) . Also intersecting at these same points are dp p /dp = and dp, n /dp = 0. Different features and behaviors 
exist between the two cases. The momentum dependent case (thick curves) has behaviors that are compressed in 
these y-p plots. The coexistence curves have a different quantitative but similar qualitative behavior between the two 
cases. The coexistence loop (dash-dotted line) is outside the other two loops and tangent to chemical instability loop 
at two points. These two points are the critical points of low and high y which are shown in FigfS] where the two 
loops touch at the peak of each loop. Comparing the two cases quantitatively, we see a compression of the results of 
the momentum dependent case (thick curves) with respect to the results of the momentum independent case (thin 
curves). The thick loops are inside of thin loops. 

FigO and Fig[H] show chemical potentials for both proton \i v and neutron /i„ and pressure P along the various 
curves of coexistence and chemical and mechanical instabilities. Curves in Fig[7| illustrate the behavior of each 
chemical potential with density on the left panel and proton fraction on the right panel. Curves in Fig[8]are pressure 
versus density on the left side and pressure versus proton fraction on the right. The separate pressure-proton fraction 
behaviors in Fig[5]were already shown in FigfSl but now these two figures contain additional plotted quantities which 
are the dp q / dp\ Vt T = curves. The chemical potential density curves in FigJTJhave features similar to those discussed 
in Figj6] Both momentum dependent and independent cases of FigJT] and Figj8] shows tangent points of the solid line 
and dash-dotted line. Also seen in these figures are the compression or shrinking of various curves for momentum 
dependent case with respect to the momentum independent case. The thick loops are inside of thin loops. FigfS] 
shows that the momentum dependent interaction leaves the point of equal concentraion nearly unchanged i.e., from 
y = 0.4106 toi/K 0.415. 

Fig[5] shows the behavior of the boundary curves of the proton and neutron chemical potentials with pressure for 
both momentum dependent and momentum independent Skyrme interactions. A comparison of the thick curves 
(momentum dependent case) and thin curves (momentum independent case) shows that the qualitative behavior 
is the same. Quantitative difference exist with the momentum independent behavior being an enlargement of the 
momentum dependent shape. The coexistence arc and the chemical instability loop meet at the cusp. The behavior 
shown in these figures also confirm that no truncation effects exist in our study. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we studied the thermodynamic properties of a two component system of hadronic matter made of 
protons and neutrons. Our analysis is based on a mean field model using a local Skryme interaction and includes 
both velocity or momentum dependent and momentum independent interactions, besides volume, symmetry and 
Coulomb effects. We have used a somewhat simplified description of the velocity dependence of the nuclear interaction. 
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FIG. 6: Figure shows the coexistence curves (dash-dotted line), chemical instability boundary curves (solid line) and mechanical 
instability boundary curves (dashed line) at T = 10 MeV. Also shown are the dfj, q /dp = curves for proton (dotted line) and 
for neutron (short dash line) at T = 10 MeV. The dash-dot-dot-dotted line is for He(p)- The thick lines are for momentum 
dependent Skyrme force and the thin lines are for momentum independent Skyrme force. The momentum dependent loops are 
inside the momentum independent loops. 

In particular we have used a density dependent effective mass approximation. Effective mass approximations are 
frequently used in physics to capture the main effects and they lead to a simpler set of equations and a corresponding 
simpler analysis. As noted we still keep Coulomb and surface terms which are present in realistic nuclear systems. It 
is the interplay of volume, surface, symmetry and Coulomb and momentum dependent terms that is studied here. In 
fact, the interplay of such terms makes nuclei a unique system for studying phase transitions, chemical and mechanical 
instability in binary systems. We then applied the basic thermodynamic relations that we developed to issues related 
to the mechanical and chemical instability of nuclei and features associated with a liquid/gas phase transition in this 
system. 

Because of the two component nature of real nuclear systems, the analysis involves a study of the behavior in proton 
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FIG. 7: Chemical potential fi p (upper panel) and \i n (lower panel) for various boundary curves at T = 10 MeV. The curves are 
same as in Fig|6] 

fraction, density and temperature (y,p,T) and also proton fraction, pressure and temperature (y,P,T). We studied 
systems with proton fraction y = ~ 1, where y = corresponds to a system of pure neutrons and y = 1 is for 
a system of pure protons. An important system with large neutron excess is a neutron star. The study of nuclear 
system with arbitrary proton/neutron ratios is also important for future RIB experiments and for medium energy 
collisions where the liquid/gas phase transition is studied experimentally. In a liquid/gas phase transition the liquid 
and gas phase have different proton fractions because of symmetry and Coulomb effects. The proton fraction in the 
liquid phase reflects a more symmetric system than the gas phase where a higher asymmetry exists. The process of 
producing a larger neutron excess in the gas phase is referred to as isospin fractionation and a review can be found 
in Ref . [l|, [2|, [J, |5|, [6|, [7| . The process is modified somewhat by the Coulomb interaction which leads to proton diffusion 
of some protons from the liquid phase back into the gas phase as discussed in Ref. [ill [ll| . 

One of the unique aspects of the nuclear systems is a velocity or momentum dependence in the two body interaction. 
Here, we also study the role of this momentum dependence first in the thermodynamic properties of the system. Then, 



20 




FIG. 9: Chemical pressure fi q versus pressure P for various boundary curves. The left pannel is for proton /%, and the right 
pannel is for neutron fi n . The curves are same as in Fig|6] 



we extend the discussion of its role to nuclear instabilities and phase transitions and make a comparison with the 
case without momentum dependence. A characteristic pattern of qualitative similarities and quantitative differences 
appear between a momentum or velocity dependent Skryme interaction and a momentum or velocity independent 
Skryme interaction. These patterns can be seen in FigQJSl and are discussed in detail in Sect lllll which we briefly 
summarizes now. 

FigQ] shows that the momentum dependence increases the pressure at a given density. Figj2]and Fig[4]show proton 
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fraction versus density and chemical potential versus density at several pressures and at a fixed temperature. The 
qualitative features arc the same between momentum dependent and momentum independent forces. However, sizeable 
quantitative differences are present between the two types of interactions. For example the solid loops (momentum 
dependent interaction) in proton fraction versus density of FigJ5] are reduced versions of the same dashed loops 
(momentum independent interaction). Similarly, the chemical instability boundaries for a momentum dependent 
Skryme interaction are found to be reduced versions of the same boundaries for momentum independent Skyrme 
interactions as can be seen from FigJS] and a comparison of the thin curves of FigJB]-[S] with the corresponding thick 
curves of these figures. Fig[5] also shows that momentum dependent terms reduce the height of the proton rich 
asymmetric loop (y > de) and leave the height of the neutron asymmetric loop (y < yE) almost unchanged while 
the lowest pressure point, which is the point of equal concentration y^, is shifted upward with the value of yE nearly 
unchanged. From FigJ5]we also see that the chemical instability loop lie on top of each other for proton and neutron as 
required by the general connection of Eq. tjoT)) . Also seen is that the chemical instability loop is inside the coexistence 
loop and tangent to it at the maxima of each loop. The largest and smallest y in the coexistence loops are shifted 
inward towards the point of equal concentration yE- Figs [6] and [7] shows that the mechanical instability loop is inside 
the chemical instability loop and tangent at the equal proton fraction yE{p) without touching it at the peak of them. 
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